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1. Introduction 

The problem of stabilizing moduli in superstring compactification has seen much re- 
cent progress. A benchmark in this progress is the work of Kachru, Kallosh, Linde and 



Trivedi |67]], which builds on many works on Calabi-Yau compactification, flux compacti- 
fication and other aspects of string compactification to propose a way to construct M = 1 
and non-supersymmetric vacua with all moduli stabilized in a controlled regime. In this 
paper, we look for specific models of the type they suggest, announce examples which we 
expect will work in all detail, and explain various senses in which these examples are less 
common than one might have supposed. 

We begin with a short review of the problem. Soon after the pioneering works on 



Calabi-Yau compactification of the heterotic string |n| , it was realized that a disadvantage 
of this construction was the presence of moduli of the Ricci-flat metrics and vector bundles 
which it uses. At least in perturbation theory, these moduli and the dilaton become 
massless scalar fields with at least gravitational strength couplings, whose presence would 
be in contradiction with experiment. 

One might think that this consideration prefers other compactifications with no mod- 
uli. However, another possibility is that the problem is just an artifact of perturbation 
theory. One can see this by postulating simple effective potentials which could plausibly 
emerge from non-perturbative effects in string theory, which have isolated minima. Any 
of these minima would be a vacuum with stabilized moduli. 

Some early examples were the models of [P3| , PT[] , and the "racetrack" models P7| , |73| , P3| 
which use non-perturbative contributions from two gauge theory sectors, with different 
dependences on the gauge coupling. Many similar constructions have been proposed over 
the years. Their essential ingredient, more than any particular feature of string theory, 
is the concept of effective potential, in which effects from various different sources are 
added into a single function which controls the vacuum structure. Once one believes 
that this is an accurate picture of the situation, then given a sufficiently rich supply of 
contributions to the potential, it becomes very plausible that the typical potential will have 
many isolated minima, so that moduli stabilization is generic. Furthermore, the value of 
the effective potential at the various minima (the "effective cosmological constant") will 
be widely distributed among both positive and negative values, possibly including the 
observed value. 

One might ask whether other general features of the problem, such as low energy 
supersymmetry, change this general expectation. As we have discussed at length elsewhere 
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PSfl , the standard expression for the potential from Af = 1 supergravity, supports this 
point of view. By postulating a superpotential and Kahler potential, one can get a wide 
range of potentials, which do not obviously favor positive or negative effective cosmological 
constant. 

One important general expectation from string theory is the following. Calabi-Yau 
compactification, and essentially any "geometric" compactification of string and M theory, 
has a "large volume limit" which approaches ten or eleven dimensional Minkowski space. In 
this limit, the four dimensional effective potential will vanish | 36fl . Recently more general 
arguments for this claim were made by Giddings pO ]. This makes it somewhat harder 
to stabilize vacua with zero or positive effective cosmological constant, as one needs the 
potential to have at least two points of inflection (the desired minimum, and a maximum or 
more generally a barrier at larger compactification radius). Still, this leads to no obvious 
problem of principle. 

To go beyond these rather general claims, one must get sufficient control over nonper- 
turbative effects to show that such potentials could arise from string/M theory, and be able 
to compute them in a wide enough range of examples to make serious predictions. This line 
of thought developed in the late 1980's, and led to a focus on nonperturbative methods in 
string theory, a field which took off with the 1989-90 work on matrix models. An equally 
important early nonperturbative result was the exact solution, including all world-sheet 



instanton corrections, of (2, 2) Calabi-Yau sigma models using mirror symmetry [20]. 

While when first formulated these techniques appeared to address only low dimen- 
sional models, the 1993-94 work of Seiberg and collaborators on nonperturbative super- 
symmetric gauge theory showed that similar results could be obtained for the effective 
superpotential in simple four dimensional models. With the 1994-95 realization [86,34,^3 
that duality and branes were key nonperturbative concepts, tremendous advances were 
made, leading to a fairly good picture of compactifications with at least 8 supercharges, 
and some success in describing compactifications with less or no supersymmetry. 

Studying such compactifications and mapping out the possibilities requires both the 
techniques to do Calabi-Yau, F theory, G2 compactification etc. and study the physics of 
each, and "classification" and "mapping" results, such as the classification of of Calabi- 
Yau's which are toric hypersurfaces, and the detailed study of their moduli spaces, for 



example \7^. These should both provide specific examples of these general recipes and 
some idea of which ones might be relevant for describing the real world. We will call upon 
a number of these results below []5T|j65| , |66| , |73|j8^ , p4 . 
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The subject of compactification with four or fewer supercharges remains somewhat 
controversial. Most works on this subject freely use field theoretic concepts such as Kaluza- 
Klein reduction and the effective potential, with the main string/M theory inputs being 
new light states and new corrections to the effective Lagrangian. At present, there is no 
real evidence that this is wrong. 

On the other hand, it is true (by definition) that a vacuum with stabilized moduli is not 
connected to the limit of arbitrarily weak coupling and arbitrarily large compactification 
volume through time-independent solutions. One must pass either through the larger 
configuration space ( "go off-shell" in an old-fashioned language) or consider time-dependent 
solutions to do this, and our lack of understanding of these subjects in string/M theory 
makes it conceivable that what look like valid solutions from an effective field theory point 
of view in fact are not valid solutions of string/M theory. This possibility has been put 
forward most forcefully by Banks and collaborators H0. While we do not find their 
present arguments convincing, the point is definitely not settled and deserves attention. 

Another issue along these lines is the question of whether there is some obstacle 
within string/M theory to realizing compactifications with positive vacuum energy. Over 
the mid-90's, convincing observational evidence was found for an accelerated expansion of 
the universe, which is most simply explained by the hypothesis that there is a non-zero 
positive "dark energy" of energy density about 0.7 times the critical value. While there 
are many possibilities for what this is, the simplest and probably easiest to reproduce 
from string theory is a cosmo logical constant, which would appear to cause the long term 
evolution of the universe to asymptote to a de Sitter geometry. From the point of view of 
field theory and the effective potential, there is no problem with getting positive vacuum 
energy. On the other hand, no such string vacuum is known, and theoretical arguments 
were even advanced that these should not exist ]6 , 45 , 54 1 . 

Meanwhile, progress in our understanding of string/M theory has led to steady 
progress in the study of moduli stabilization. An important recent work in this direc- 
tion is that of Kachru, Kallosh, Linde and Trivedi [67], who argue that metastable de 
Sitter vacua with all moduli stabilized can be constructed in type lb string theory, at least 
in an effective potential framework. 

The KKLT work has two aspects. In the first, they propose a recipe for such a 
construction, which we will discuss in detail shortly. In the second, they sidestep the claim 
that de Sitter space could not be realized in string theory, by suggesting that our universe 
need not asymptote to de Sitter space. Rather, our vacuum might be metastable, in a 
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positive local minimum of the effective potential. As we dicussed, since the potential goes 
to zero at large compact ificat ion volume, any such vacuum is potentially metastable to 
decay to this decompactification limit. Using results of Coleman and de Luccia, Hawking 
and Moss, and others on tunneling in semiclassical quantum gravity, KKLT argued that for 
almost any reasonable potential, one would find an extremely small but non-zero decay rate, 
slow enough to pose no cosmo logical constraint, but fast enough to evade the paradoxes 
which had been suggested might forbid eternal de Sitter space as a vacuum. These ideas 
are reviewed in ||] (see also the recent H51fl). 

This leaves the more technical problem of showing that such metastable de Sitter vacua 
actually do exist in string theory compactification. KKLT suggested a recipe by which this 
could be done, following previous work of many authors on moduli stabilization by fluxes (a 
very incomplete list includes ||g^g0^^7lg7|,§§§g§§,§g|§§|) and especially Giddings, 
Kachru and Polchinski KJ who developed a class of lb flux compactifications which sta- 
bilize the dilaton and complex structure moduli. Some advantages of this construction are 
that it allows fixing the dilaton to weak coupling by making an appropriate choice of flux, 
and that cancellations can bring the flux contribution to the potential below the string 
scale. Most importantly, it is based on the relatively well understood theory of Calabi-Yau 
moduli spaces, so explicit computations can be done. 

Having this construction in hand, the remaining problem is to stabilize the Kahler 
moduli and obtain a small positive cosmo logical constant. Granting the effective poten- 
tial framework, one proceeds as in the early works on moduli stabilization, to look for 
nonperturbative contributions to the superpotential which could stabilize the remaining 
moduli. As it happens, in lb compactification nonperturbative effects from gauge theories 
on branes, and from D-instantons, all depend on Kahler moduli, and thus one would expect 
that generically this stabilization would not be a problem - once one had computed these 
constributions to the effective potential sufficiently well, they would stabilize all moduli. 

Finally, to obtain a positive cosmological constant, one has several choices. Con- 
ceptually the most straightforward is to look for supersymmetry breaking vacua of the 
previously computed effective potential, which again on grounds of genericity should exist. 
Alternatively, one can look for supersymmetric AdS vacua which stabilize all moduli, and 
then add other supersymmetry breaking effects which lift the vacuum energy. This idea 
has the advantage that the resulting potential can naturally have the barrier required for 
metastability of a de Sitter vacuum. KKLT suggested to add anti D3-branes, and one can 



imagine many other possibilities such as D terms coming from brane sectors [p9|Jlq , |12 
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The upshot is that such effects, if they exist, could well combine to stabilize all moduli 
at a metastable de Sitter minimum. There are two obvious gaps in this argument. The 
deeper one is the question of whether the effective potential analysis is valid. Perhaps 
the best way to address this is to consider a specific example, and try to justify all of its 
ingredients in ten dimensional string theory. 

To do this, one must address the first gap, which is that no concrete model of this 
type has yet been put forward. In this work, we will fill this gap, by examining models of 
the general class discussed by KKLT, to see if this idea can actually be implemented in 
string theory. We will need to call upon many of the ideas discussed above, especially the 
classification of Calabi-Yau three-folds and four-folds, and the analysis of nonperturbative 
effects in F theory by Grassi [57]. 

While we will exhibit models in which nonperturbative effects lift all Kahler moduli, 
we will also argue that, at least if we rely on instanton effects, such models are not generic. 
There is a fairly simple reason to expect this. In the language of orientifold compactification 
with branes, it is that these world-volume gauge theories usually have too much matter to 



generate superpotentials, as seen in works such as [[15), and as we will see in a large class 
of models below. This is because the cycles on which the branes wrap can be deformed, 
or because the branes carry bundles with moduli. In either case, the gauge theory has 
massless adjoint matter, which eliminates the superpotential. The simplest case is the 0- 
cycle or point, which clearly always has moduli. Of the two-cycles, only S 2, s can be rigid; 
higher genus Riemann surfaces in Calabi-Yau almost always come with moduli. Bundles 
on surfaces (and the entire CY) normally have moduli as well. 

On the other hand, in a sizable minority of models, the instanton generated superpo- 
tentials are sufficiently generic to stabilize all Kahler moduli. The examples that we can 
check systematically are toric fourfolds X elliptically fibered over a three dimensional base 
B which is Fano, and P 1 -bundles over a toric surface. Much of what we say should hold 
more generally. 

There are no very simple models in this class. In particular, one can argue that no 
model with one Kahler modulus can work.0 

Of the models which work, we have considered three of the simpler ones in some 



detail. The fourfolds are elliptic fibrations over the Fano threefolds T\\ and JF 18 [S3, 84, 79 



1 This conclusion was reached independently by D. Robbins and S. Sethi [|8|. Also, relevant 
work of V. Balasubramanian and P. Berglund is mentioned in [H . 
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Their Euler characteristics are x = 16848 and x = 13248 respectively. Some of our 
considerations will require that we work with their orientifold limits. These are Calabi- 
Yau threefolds that can be realized as hypersurfaces in toric varieties and their Hodge 
numbers are h ' = 3, h 2,1 = 111 and h > = 5, h 2,1 = 89 respectively. Finally, we consider 
an elliptic fibered CY over P 2 , a much studied CY El with h 1 ' 1 = 2 and /i 2 ' 1 = 272. 



2. Review of KKLT construction 



We start with lb superstring theory compactified on a Calabi-Yau orientifold, specified 
by a choice of CY threefold Z and a holomorphic involution fl on Z whose fixed locus 
consists of points and surfaces (03 and 07 planes). For example, we might take Z to be 
the hypersurface in CP 4 defined by a quintic polynomial f(z) = 0, and Cl to be the map 
z\ — > —Zi, Zi — > Zi for % > 1. 

We then introduce D-branes in such a way as to cancel the RR tadpoles produced by 
the orientifold planes, as in [[52] . For a recent discussion of this problem, see |16| . 



The number of complex structure moduli of Z is h 2,l (Z), and the number of complex- 
ified Kahler moduli is h l,l {Z). We are going to refer to these as "shape" and "size" moduli 
respectively, mostly because the words "complex" and "Kahler" have so many different 
roles in the discussion that the usual terminology is cumbersome (one gets tired of saying 
"Kahler metric on the Kahler moduli space," etc.) 

A more general starting point would be F theory ||92| . We will use this language 
to construct models below, but rapidly move to the Po orientifold limit in our examples, 
since there are many unanswered physical questions in either picture, which we will need 
to appeal to the underlying definitions in string theory to resolve. 

The basic relation is as follows. An F theory compactification is defined by a choice 
of Calabi-Yau fourfold X with an elliptic fibration structure over a threefold B, 



->■ X 



(2.1) 



B. 

Physically, we compactify Po theory on B, introducing 7-branes at the singularities of the 
fibration it, so that the resulting dilaton-axion field at a point p G B corresponds to the 
complex structure modulus r of the fiber 7r _1 (p). The orientifold limit [|87| is then the 
special case in which all of these singularities are singularities (an 07-plane and four 
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coincident D7's). In this case, Z is a double cover of B branched at the singularities, and 
Q exchanges the two sheets. Let us denote the part of the cohomology of Z even or odd 
under Cl with subscripts ±. 

Classically and in the absence of fluxes, this compact ificat ion has a moduli space of 
Af = 1 supersymmetric vacua, parametrized by hi shape moduli, the dilaton-axion, h 1 ^ 1 
size moduli complexified by the RR 4-form potential C4, and hi moduli from the 2-form 



potentials B2 and C2 [|T5| J 

We will avoid having to discuss the 2-form moduli by restricting attention to models 
with hi = 0. More generally, while we do not know a mechanism which would stabilize 
them at large volume, they might be stabilized by world-sheet and Dl-instantons, or by 
couplings to brane world- volumes. 

In general, there are also open string moduli corresponding to positions of D7-branes 
and D3-branes, and moduli of bundles on D7-branes. We will ignore these through most 
of the discussion and return to them in the conclusions. One excuse for this is that these 
moduli spaces are all expected to be compact. This is intuitively clear for positions of D3 
branes, since their moduli space is Z itself, which is compact. This should also be true for 
moduli spaces of bundles except for small instanton limits, but these just correspond to 
other brane configurations. Thus, there is no analog of the decompactification or runaway 
problems for these moduli, and generic corrections to the potential should fix them. 

2.1. Fixing moduli using fluxes 

We begin by trying to fix the shape moduli and the dilaton-axion. In lb language, this 
is done by turning on the NS and RR three-form field strengths, and respectively. 
As discussed in many references (e.g. [f49|,|2|1) the equations of motion will force these to 
be harmonic forms, so they are determined by their cohomology classes in H 3 (Z, IR). 
There is a Dirac-type quantization condition which normally forces these classes be integer 
quantized in units of the string scale. Setting this unit to one, they live in H 3 (Z, 

The choice of flux is constrained by the tadpole condition on the RR four-form po- 
tential, 

L=\J G 4 AG, = ^rj a0 N NS = ^-N m , (2.2) 

where Nds is the number of D3-branes minus the number of anti-D3 branes. 

The ten-dimensional lb supergravity analysis of unbroken Af = 1 supersymmetry in 



this context is done in |p6[ ; the most important condition is that 

G (3) = p(3) _ tH (3) 
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must be imaginary self-dual, G = i * G. One can restate this as the condition that 

G (3) G #(O,3)( Z)C ) 0iJ (2^ ZjC )_ 

Setting the remaining parts of the cohomology to zero amounts to h 2,1 + 1 conditions, 
which is precisely enough to fix all shape moduli and the dilaton. Thus, we might expect 
generic choices of flux to do this. 

There is an additional condition for this to be true with four noncompact Minkowski 
dimensions: the _ff"(°' 3 ) component of G^ must be zero. This is one more condition than 
the number of moduli, so this is non-generic. Otherwise, we get supersymmetric AdS 
vacua. 

A nice physical summary of these results, which will be essential later on, is that these 
conditions on G^ follow from solving the supersymmetry conditions in an Af = 1 effective 
supergravity theory. Its configuration space is the combined Calabi-Yau and dilaton-axion 
moduli space; in other words the product of the shape moduli space M. C (Z), the size 
moduli space A4k(Z), and the upper half plane 7i, with the Kahler potential 

K= -loglmr-log / OAO-21ogF (2.3) 

where r is the dilaton-axion, V is the volume of the Calabi-Yau as a function of the size 
moduli (more on this below), and O is the holomorphic three- form on Z. This is just 
the Kahler potential of the compactification with zero flux and M = 2 supersymmetry, 
computed in the a', g s — > limit. 

As superpotential, we take the Gukov-Vafa-Witten superpotential 

Wqvw = [ ft AG (3) . (2.4) 



One can see [f|9],|2| that the conditions 

= DiW = diW + {diK)W (2.5) 

for the shape moduli set the H^ 1 ' 2 ' part of G^ to zero, and for the dilaton-axion sets the 
_£/"( 3 '°) part to zero. 

For typical values of the flux, this will fix r and the shape moduli z at a mass scale 
mo ~ a'/yV, where V is the volume of Z. The typical value of e^lW] 2 at the minimum 
is 1 in string units, but this can be made smaller by tuning the fluxes. One expects its 
smallest value to be ~ l/A/"„ ac where M vac is the number of flux vacua on the complex 
structure moduli space, as we discuss below and in |2^ . 



We will discuss the technicalities of finding choices of flux which stabilize these moduli 
in a desired region of moduli space in section 4. It will turn out that this is often, but not 
always possible. 
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2.2. Fixing size moduli 

Since ( |2.4j ) is independent of the size moduli, the conditions DiWcvw = for the 
size moduli amount to either diJC = for all size moduli, or else W = 0. It is easy to 
see (as we discuss below) that in the large volume limit, one cannot solve diJC = for all 
moduli. One can find Minkowski supersymmetric solutions with Wgvw = 0. One can also 
find Minkowski but non-supersymmetric vacua with Wgvw 7^ 0, because of the no-scale 
structure of ( |2.3|) . 

In either case, while we might find other effects which lead to a positive vacuum 
energy, this energy will typically decrease with increasing compactification volume, and the 
resulting solutions will be unstable. We need a more complicated effective potential to fix 
this problem, which might be obtained by incorporating stringy and quantum corrections 
to W and K. 

Unfortunately doing this remains a hard problem, both to get results in examples, 
and even to define what one means by the nonperturbative effective potential in general. 
In particular, there are essentially no results on nonperturbative corrections to the Kahler 
potential. 

To address this, one can seek examples where the instanton expansion appears to be 
valid, so that the leading corrections will dominate. This was the primary goal in KKLT's 
discussion and we will address it below. However we believe that there is no question of 
principle which requires restricting attention to these examples; rather we do this because 
of our limited technical control over the theory at present. 

KKLT suggested to stabilize the size moduli using quantum corrections to the su- 
perpotential, as they are known to be present and depend on the size moduli in many 
examples, and are controlled at large volume and weak coupling. 

For example, one can consider a case in which D7-branes wrapped on a cycle E have 
pure SU(N) Yang-Mills theory as their world- volume theory. The gauge coupling in this 
theory is \ = V(E) where V(E) is the volume of the cycle E in string units, and thus 
nonperturbative effects in this theory lead to a superpotential 

W eff = Be- y ^/ N (2.6) 

where B is a (presumably order 1) quantity determined by threshold effects, etc. and 
which may also depend on size and shape moduli. 
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If we consider a vacuum in which the flux contribution to the superpotential ( |2.4| ) 
takes a value Wo independent of size moduli, and grant that this is the leading correction, 
the supersymmetry condition for the size modulus becomes (let p = iV), 

= D P W = - — e lp l N - — (Wo + Be^ N ) . (2.7) 
N Imp V / 

For small Wo, this will normally have a unique solution at large imaginary p, roughly given 
by 

, 3NW 
p ~ -i log - . 

To see this without detailed calculation, note that the problem of solving D p W = with 
W 7^ is the same as the problem of finding critical points of the function / = e^lW] 2 . 
This function is positive and goes to zero at very large Imp. If we contrive W to have a 
zero at large Imp, then it is easy to see that as we increase Imp from this value at fixed 
Rep, we are always at a critical point in Rep. Furthermore, as we increase Imp, / starts 
out increasing, so for it to later decrease it must pass through a critical point in Imp. 

This solution will have W ~ Wo and thus we will end up with a supersymmetric AdS 
solution with small negative cosmo logical constant, controlled by the previous step where 
we chose fluxes to get small Wo- 

Another source of similar non-perturbative contributions to the superpotential is D3- 
brane instantons wrapped on surfaces in Z ||94|| . Their dependence on the size moduli is 
the same and thus this level of the discussion works the same way. As we discuss below, in 
the case that D7-branes wrap a cycle E, these are equivalent to gauge theory instantons, 
but are present more generally. Thus in the detailed discussion below, we consider the 
non-perturbative corrections as generated by D3-brane instantons. 



Some further general observations on this problem appear in [17 



2.3. Breaking supersymmetry 

We will have enough trouble trying to realize the first two steps, but given this, to go 
on to the third step we would need some approximate expression for the warp factor on 
Z , so that we could argue that the increase in vacuum energy produced by adding an anti 
D3-brane could also be made small. The simplest way to imagine this working is to fix the 
shape moduli near a conifold singularity, and appeal to the Klebanov-Strassler solution 
||72j| as an approximate description of the metric in this region. 

While it seems to us that this should generally work, of course there is room for more 
subtleties, as we will briefly discuss below. 
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3. Details of size moduli stabilization 

We seek the conditions on Z which will lead to a superpotential which stabilize all 
size moduli significantly above the string scale, so that ot corrections can be ignored. It 
will emerge from considerations below that there are no models with h ' = 1, and thus 
we consider the multi-moduli case from the start. We will usually set £ s = lT{\fa' to 1 in 
the following. 

Let Di be a basis of divisors on Z (essentially, these are classes in H±{Z, Z)). The 
main data we will need about Z are the triple intersection numbers of the divisor basis, 
Dij k = Di-Dj-D k . Using Poincare duality, we also let Di denote the corresponding class in 
H 2 (Z, Z). For CY's which are hypersurfaces in toric varieties, there are efficient methods 
for computing these intersection numbers and the other data we are about to discuss. We 
discuss some relevant examples in section 5 and appendix B. 

We define the size moduli by writing the Kahler class J of Z as 

J = J2t i D i , (3.1) 

i 

defining a set of real coordinates t % on the space of Kahler classes. 

The classes which actually correspond to Kahler metrics are those lying in the Kahler 
cone, defined by 

0< [ J (3.2) 
Jc 

for all holomorphic curves C. The classes of such curves are called "effective classes" and 
form a cone, the Mori cone. 

In fact, to be able to ignore a' corrections, we want all such areas of curves to be at 
least 0(1) in string units. 

The natural holomorphic coordinates which appear in lb orientifold compactification 
are not the t^, but are instead the complexified volumes of divisors, which we denote r^, 

n = [ ljAJ-iC 4 . (3.3) 

JDi ^ 

We denote their real parts as 

Vi = Re n. 

The Kahler potential on this configuration space, neglecting a' and g s corrections, is 

K K = -2 In V(t,t) (3.4) 
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where V is the volume of Z, defined by 



V = - J 3 = -D ijk tHH k . (3.5) 
6 6 3 K ' 

We can write this (implicitly) as a function of r» using 

n = w = \ DiJ2 = \ D ^ ktHK ( 3 - 6 ) 

This change of coodinates is essentially a Legendre transform on V, as discussed in detail 
in 
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3.1. Instanton superpotentials 

The best studied corrections to the superpotential in F theory are produced by D3- 
brane instantons wrapped around divisors. These take the form 

W np = Y,bne- 2M (3.7) 

n 

where the b^ are one- loop determinants on the divisors = n l Di, depending on the 
complex structure moduli. At large radius, 2-KTi ^> 1 and the mass scale of this contribution 
to W is exponentially suppressed compared to mo, so to a good approximation we can 
consider the complex structure moduli to be fixed by the tree level flux superpotential Wq, 
and use W = Wo + W np with constant Wq and b^ as effective superpotential for the Kahler 
moduli. 

The most important thing to know about the determinants is whether they are 



non-vanishing. This was studied in |53| by using the relation between F theory and M 
theory compactification on the four-fold. In M theory, such a correction comes from a 
fivebrane instanton wrapped on the lift Vfi of the divisor to the fourfold. The F theory 
limit is the limit in which the area of the fiber goes to zero, and a divisor V will contribute 
in this limit if it is vertical, meaning that 7r(V) is a proper subset of B. 

In this formulation, b^ includes the determinant of a Dirac operator on V, which can 
have zero modes. In fact the zero modes turn out to be equivalent to holomorphic p-forms 
on Vfj tensored with a three-dimensional spinor. Let the number of such p- forms be h 0,p . 

An instanton which contributes to W np must have two fermion zero modes (corre- 
sponding to the two supersymmetries broken by the instanton). Since h°'° = 1, this will 
be true if h 0,J = for all j = 1,2,3. On the other hand, it is possible for other brane 
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couplings (in particular, couplings to the NS and RR fluxes we used to stabilize other 
moduli) to lift zero modes in pairs. Thus it is possible that b^ ^ more generally, but in 
any case the necessary criterion for this is that is a divisor of arithmetic genus 

x {O v ) = J2(-l) j h°> j = I- 

j 

In this more general case, it is also possible for variations in complex structure to bring 
down or lift zero modes in pairs, so b^ could have zeroes at special points in complex 
structure moduli space, possibly depending on the choice of flux. 

Unfortunately not much more is known about the determinants b^. Note however 
that they must be sections of the same line bundle as the holomorphic 4- form O4 on J; 
otherwise the combined superpotential wouldn't make sense. Let us fix the overall scale 
of O4 (and thus 6^) by putting J x Q4 A O4 = 1 at the given point (r, z). In this natural 
normalization, it is reasonable to assume the are generically not much smaller or bigger 
than 1. 

The upshot is that divisors whose lift to X have = for all j = 1,2,3 will 
contribute to the superpotential, whether or not flux is turned on and regardless of the 
values of the other moduli. The most general class of divisors which can contribute are 
the divisors of arithmetic genus 1; these contributions might vanish for special values of 
flux and mo dulii 

3.2. Gauge theory superpotentials 

Superpotentials produced by non-perturbative effects in gauge theories living on D7- 
branes can also be understood as coming from divisors of arithmetic genus 1. 

Roughly, if one wraps a D7 about a divisor E, the zero size limit of the Yang-Mills 
instanton in the resulting world- volume gauge theory is just a D3-instanton wrapped on E. 
A detailed discussion is somewhat more complicated. For example, in pure SU (N c ) super 
Yang-Mills theory, the instanton has 2N C fermion zero modes, which for the corresponding 



One might worry about possible modifications to this analysis if V is spin c and not spin, 



for reasons discussed in [46|. These can modify the coefficient b and might even introduce dilaton 
dependence |^2| , but do not affect the statements made here or our later conclusions. We thank 
G. Moore and S. Sethi for communications on this point. 

We also note that the analysis in [94J used a U(l) symmetry of the normal bundle to the 



divisor. If this were broken by the fluxes, the arithmetic genus 1 condition might be relaxed [70]. 
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D3-instanton arise from 3 — 7 open strings. Thus it does not contribute directly to the 
superpotential. Rather, its effects force gaugino condensation which leads to a vacuum 
energy which can be reproduced by a superpotential. This accounts for the 1/N in the 
exponent of (P75|) . 

In F-theory, nonabelian gauge symmetries arise from singularities of the elliptic fibra- 
tion. For example, an A^v-i singularity fibered over a surface S in the base corresponds to 
an SU(N) gauge theory living on S. If S has h ' 1 = h 0,2 = and the singularity does not 
change over S, there will be no additional matter and we are in the pure SU(N) case just 



discussed. In [71] , it was shown that by compactifying on a circle and going to the dual M- 
theory picture, the three dimensional nonperturbative superpotential can be computed as 
coming from M5 brane instantons wrapping the exceptional divisors obtained by resolving 
the An_i singularity. Again, these divisors are of arithmetic genus 1. The correct four 
dimensional result is then obtained by extremizing the three dimensional superpotential 
and taking the decompactification limit. 

A similar geometrical analysis can be done when in addition fundamental matter is 
present [^2j. The situation is more subtle here, but the upshot is that the gauge theory 
superpotentials can again be derived from M5 brane instantons wrapped on divisors of 
arithmetic genus 1. 

Thus, in the known examples, a necessary condition to generate a nonperturbative 
superpotential is the presence of a divisor of arithmetic genus 1 in the fourfold.! 

A more intuitive but only partial argument for the role of the arithmetic genus in the 
brane gauge theory is that the cohomology groups H°' P (V) which entered our previous 
discussion also control the matter content on the D7-brane gauge theories. In particular, 
non-zero h 0,1 and h 0,3 will lead to massless adjoint matter, which tends to eliminate the 
superpotential (as in Af = 2 gauge theory). This conclusion is not strict as gauge theory 
with a non-trivial world- volume superpotential for the adjoint matter (say Tr <p n ) can 
generate a superpotential, but this possibility comes along with the possibility of deforming 
the superpotential to give mass to the adjoint matter and thus is accounted for by the 
possibility we discussed earlier of lifting the zero modes in pairs (this can be made more 
precise by relating the superpotential to obstruction theory and h 0,2 as in the D-brane 
literature 



3 Actually, in [94], a more general possibility is raised: divisors of arithmetic genus x(-D) > 1 
(but presumably never X < 0) might contribute to the superpotential, through strong infrared 
dynamics or "fractional instanton" effects. We thank S. Kachru for reminding us of this point. 
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3.3. Complete sets of divisors 



Divisors of arithmetic genus 1 tend to be rare [|94] , |57H . Furthermore, in order to stabilize 
the compactification at strictly positive radii, a sufficient number of distinct divisors must 
contribute. This can be seen as follows. 

We are looking for solutions to D T .W = d T W + {d T JC)W = 0. Now 

dJC K = -±J 2 dJ = -^JdJ 2 = (3.8) 



hence d Ti IC = — ^7, and 



= D T W = Y J bne- 2 ™- 7 {-n l ) - —W. (3.9) 



By multiplying with Di and summing, it follows that 



J^bne-^Drt + ^-J = 0. (3.10) 



We want vacua with W 7^ 0. This is generic, and if the contributing divisors are linearly 
independent, as is often the case in examples, it follows directly from this equation. 

In this case, the Kahler class J is a linear combination of the contributing divisors. 
Thus, a necessary condition for a solution with positive radii to exist is that a linear 
combination 

7? = r-D- 

of the contributing divisors exists such that 1Z lies within the Kahler cone (|3.2| ). We will 
call such a set of divisors "complete." 

Note that a set of divisors cannot be complete if there exists a curve C which does not 
intersect any divisor in that set, because then all 1Z will lie outside the Kahler cone, or at 
best on its boundary. On the other hand a full basis of divisors is automatically complete. 
Since it is easy to verify if a set of divisors is a basis, we will focus mainly on examples in 
which a full basis in H^Z) of divisors contributes to the superpotential. 
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3. 4- Parameters 

At large volume, the exponentials in ( |3.10| ) are very small, so for a solution to exist, W 
and therefore Wq has to be very small, of the order of the dominant exponentials e _27rn T . 
Again because J in ( |3.10| ) has to lie well inside the Kahler cone, the set of contributing 



divisors which have e _27I " n r ^ Wq has to be complete in the sense introduced above, 
since the terms for which e~ 27Tn ' T <C Wq can effectively be neglected. 

For a suitable choice of basis of divisors, the t % in (|3.1| ) correspond to curve areas, so a 
rough estimate for the divisor volumes is V* ~ kA 2 , where A is the average curve area and 
k some constant, which increases with the number of Kahler moduli due to the increase 
in terms in ( |3.6|) (in examples k is at least proportional to the number of Kahler moduli). 
So an estimate for the average curve size at a critical point of the superpotential is 

We see from this estimate that for many Kahler moduli (so k is large), Wq has to be 
exceedingly small to have A > 1. 

Moreover for many moduli the probability that some curves are significantly smaller 
than the average increases, making it even more difficult to fix the moduli in a controlled 
regime. This motivates us to look for models with few Kahler moduli. 

3.5. Geometric considerations 

There are some general geometric constraints on the existence of divisors of arithmetic 
genus 1. First, we can assume without loss of generality that there exists a Weierstrass 
model 7Tq : W — > B and that \i : X — > W is the resolution of the Weierstrass model. 
We recall that in F theory only the vertical divisors of arithmetic genus one contribute to 
the superpotential. As explained in ||57|| , such divisors fall in two classes: they are either 



components of the singular fibers or they are of the form V = iz*(D) for D a smooth divisor 
on B. 

To construct elliptic fourfolds which have divisors of the first type, one can take B to 
be a P 1 -bundle over a surface B' and enforce an ADE singularity of the Weierstrass model 
along B' ', as in ||71|| . The exceptional divisors of \i : X — > W will have arithmetic genus one. 
In lb string theory terms, this corresponds to realizing pure Yang-Mills theory with ADE 
gauge group on 7-branes (D7-branes in the A and D cases, and for the E case as in fl26j) 
and wrapping these branes on B' . The arithmetic genus one condition then corresponds 
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to the condition h 0,1 (B') = h 0,2 (B') = for the resulting four dimensional gauge theory 
not to have matter. Because of the association of these divisors to brane gauge theory we 
refer to them as "gauge-type divisors." 

The other class of divisor, pullbacks of smooth divisors on B, can give rise to D3- 
brane instanton corrections which need not have a gauge theory interpretation, and so 
we refer to these as "instanton-type divisors." In p^ , [78| , |73|j79|1 it was proposed, based on 
the study of examples, that such divisors are always "exceptional" in the sense that there 
exist birational transformations of Calabi-Yau fourfolds which contract these divisors. This 
was shown to be true whenever B is Fano in fl57|. Moreover, when B is Fano or toric, the 



number of the divisors contributing to the superpotential is finite. This is because they are 
the exceptional divisors associated to the contraction of one of the generators of the Mori 
cone, which is polyhedral for B Fano or toric.@ These divisors have a negative intersection 
with the corresponding generator of the Mori cone and are thus "non-nef," as first pointed 
out in 



One can now see that there are no models with h ' = 1, because the negative inter- 
section condition in this case reduces to 

E-D < 

where D = 7r(V) and £ is an effective curve.@ The positivity condition ( |3.2[ ) then forces 
the Kahler form to be J = —tD with t > 0. The conditions V = J 3 /6 > then forces 
Vol(D) = D ■ J 2 /2 < 0, so such an instanton correction cannot exist. In words, the 
"exceptional" nature of these divisors means that in each of the instanton amplitudes 
( p.7| ), at least one of the coefficients in the action n ■ f must be negative (in some basis). 
Given several Kahler moduli and divisors, this need not be a problem, and we will find in 
examples that it is not, but it does preclude h ' = 1. 

A more mathematical argument for this point is that if arithmetic genus one divisors 
are indeed always exceptional, then there are no models with h > = 1, because contract- 
ing a divisor will decrease h ' (from the "Contraction Theorem" cited in [57]), which is 
impossible in this case. 



It is also shown in [|57| that for B Fano, divisors of arithmetic genus one in fact have 



h°' p = for p > 0, so the necessary condition for an instanton correction is in fact sufficient. 



There are other examples in which the number of contributing divisors is infinite, such as 
the example of 



5 This inequality can be violated for gauge-type divisors, but we believe all models containing 
these have h 1 ' 1 > 1. This point is discussed further in fl8q|. 
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4. Details of flux stabilization 

In the next section, we will find that models which will stabilize all Kahler moduli 
have many complex structure moduli, n ~ 100. We will then need to find choices of 
flux which stabilize the other moduli and lead to small Wq. While one can search for 



solutions numerically, as done in previous work [68,53,81], computation time becomes an 
important issue. In particular the most straightforward approach of picking arbitrary flux 
vectors, trying to find solutions to DWq = 0, and hoping the solutions satisfy the desired 
properties, especially that of having small but non-zero Wq, becomes infeasible. One can 
simplify the problem somewhat by imposing discrete symmetries, as we discuss below, but 
this of course ignores most of the possible vacua. 

A simpler goal is to obtain information about the existence and number of flux vacua 



using the indirect approach of [p| ^9|j44| . We want to know how many supersymmetric vacua 
we can expect with L < L*, | VKo | 2 < A*, the dilaton r within a region H C H and the 
complex structure z within M. c C Ai c , in the limit of very small A*. By approximating 
the sum over flux lattice points by a volume integral in continuous flux space, this was 
computed in |2^] to be 



N vac = ( 27rL *) &3 ^L vo l(7i) / d 2n z det s po{z) (4.1) 

where 63 = 2n + 2, g is the Weil-Petersson metric on A4 C , and po a certain density function 
on M. c computed from local geometric data. The detailed expression for p and a discussion 
of its evaluation for large n is given in appendix A. A useful estimate is the index density 

detgpo ~ ^detiR + u) (4.2) 

with R the curvature form and 10 the Kahler form on Af c . In particular po 7^ 0, so the 
distribution of vacua with [Wo] 2 < A* is uniform in A* around A* ~ 0. 

Performing the integral Q4.1| ) over a region of moduli space provides an estimate for the 
number of quantized flux vacua in that region. While the estimate only becomes precise 
in the limit of large L*, numerical experiments suggest it is fairly good for L* > 63. In this 



case, one expects a subregion of radius r > ^Jb^jL^ with an expected number of vacua 
N vac ^> 1 to contain flux vacua [|9 . 



Once we know flux vacua exist in some region, there are better ways to find explicit 
flux vacua. We have developed a method which begins by fixing a rational point at large 
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complex structure, finding the lattice of fluxes solving DW = and finding short lattice 



vectors using advanced algorithms [^5J . One can then move in by systematically correcting 
the point on moduli space to take into account the corrections from this limit. Since these 
corrections are small, this often produces vacua with small Wq, and we will cite some 
results obtained this way below. 

4-1. Metastability 

Although we will not study the question of whether one can break supersymmetry 
by antibranes or D terms in any detail, it leads to another constraint on the flux vacua 
which we will study: namely, we must insist that the potential is actually minimized at 
the candidate vacuum. This was not required for consistency of a supersymmetric AdS 
vacuum [14], and is trivial for a no-scale nonsupersymmetric vacuum, but this condition 



becomes non-trivial after size modulus stabilization. 



As discussed in |29| and many other references, the mass matrix for a vacuum satisfying 
DW = 0, and for fields which do not participate in the D-type supersymmetry breaking, 
is 



m 2 = H 2 - 3\W\H, 
where we defined W = e K / 2 \W\, and the matrix 

H = 2d 2 \Wl 

expressed in an orthonormal frame. Thus, positive eigenvalues of H which are less than 
S\W\ will lead to tachyons. In the KKLT construction, Wq is assumed small, and we show 
below that this implies that \W\ is small, so that this need not be a stringent condition. 
We again take the superpotential to be a sum 



W = W (z 1 T) + W p (p). 



(4.3) 



The matrix H is given by 



H = \W\ ■ 1 + 



with 





(° 


s 





T \ 


1 


s 





T 





\w\ 





T t 





u 




\ft 





u 


0/ 



S = WD.DjW, U = WD a D a W, T = WDjDnW 
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where i, j are ortho-normal frame indices for the complex structure moduli and the dilaton, 
and a, (3 for the Kahler moduli.! If S ^> U, T, as is generically the case in the KKLT 
construction, the effect of T on the eigenvalues of H will be subleading. This expresses 
that when there is a large scale separation of complex and Kahler moduli masses, they can 
be treated separately. Moreover in general, if tachyons are found when T is set to zero 
(i.e. purely in complex structure or Kahler directions), there will also be tachyons in the 
full problem. Thus it is a good and useful approximation to put T = 0, and we will do so 
in the following. 

The condition on the complex structure moduli is then the same as in the discussion 



of [^9[ in which Kahler moduli were ignored, except that W is shifted to include W p . This 
sets the overall scale of the condition on D 2 Wfi ux but does not enter in the details. For 
the simplest model superpotential 

W p = be tp/N , 

we have 

11 =11 -3 — ] Wn 



Wn 

l + 2Imp/3NJ 

at D P W = 0, so for typical p and N this decreases W and makes tachyons less likely, but 
not dramatically so. 

For generic mass matrices, in a model with n complex structure moduli, some of these 
would be tachyonic in roughly a fraction 3n|VF| of cases, which is small for small \W\. 
Going on to consider the Kahler moduli, for one Kahler modulus, there are no tachyons 

if 

g PP\DlW\ > 2\W\ 



which is 



4 (Imp)' 



d 2 D W p — - d p Wp - ^ , W 

p p 2Tmp p p 4(Imp) 2 



> 2IWI. 



6 It may be counterintuitive that the mixing T is present as naively both K and W are the 
sum of two independent functions, so the two sectors appear to decouple. One way to see why this 
is naive is to note that this statement is not invariant under Kahler- Weyl transformations, even 
those which take the factorized form W — > fo(z,r)f p (p)W. A more mathematical way to say this 
is that the decomposition K = K(z, r) + K(p) implies that W is a section of a tensor product 
line bundle, and in writing ( }4,3| ) one is implicitly choosing reference sections of the two bundles. 
The mixing matrix T is then — \W\ 2 times the tensor product of the covariant derivatives of these 
reference sections. 



20 



For the model superpotential this becomes 



1 + 



Imp 



> 1 



which is always satisfied. This might fail in a multi-modulus model, though the examples 
we study below were found numerically to be tachyon-free as well. 

Thus, in general one does not expect moduli to become tachyonic after supersymmetry 
breaking. However special structure in the mass matrix might change this conclusion, and 
the main point of this discussion is that (in the approximation that we ignore the one-loop 
determinants in W p ) the mass matrix for complex structure moduli can be analyzed in 
the simpler model in which Kahler moduli are simply left out of consideration, since their 
effect is just to renormalize Wo- This justifies the analysis of [p|,[29[] in which this was done. 

In [p9[| , the one parameter models were studied in great detail, and it was found that 
tachyons are generic in some regimes, for example near conifold points. This is potentially 
important as we might want to work near a conifold point to obtain a small scale or small 
supersymmetry breaking. The situation for multi-modulus models appears to depend on 
details of the specific model, and this might or might not be a problem. 



5. Search for models 

The upshot of the previous section is that a model which stabilizes Kahler moduli must 
be based on a fourfold X such that the divisors on Z whose pullbacks have arithmetic genus 
one form a "complete set" as discussed in subsection 3.3. The simplest way this can happen 
is if such divisors form a basis of H^(Z). We now show this is a rather strong requirement 
and that not very many models satisfy it. 

There are two large classes of examples we consider: Fano threefolds and P 1 bundles 
over toric surfaces. 

First, since there is a classification theory of Fano threefolds, all such examples can 
be listed. In particular, |57] lists all the Fano threefolds together with the divisors that 



give rise to divisors of arithmetic genus one in the associated elliptically fibered fourfold. 



Out of all the toric threefolds in the tables in [57] , the models with a basis of divisors 
of arithmetic genus 1 are B = jF n , JF 12 , JF 14 , JF 15 , JF 16 and J-'is- The model with B = Tyj 
comes close, but while it has enough fourfold divisors with arithmetic genus one, the 
corresponding divisors in the base do not generate the Picard group of Tyj. Thus, out 
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of the 18 toric Fano manifolds in the list, only six work. We note that there are also 74 
Fano threefolds that are not toric. Out of these, only 23 have enough fourfold divisors of 
arithmetic genus one. 

We consider JF 18 and T\\ in more detail below, because their orientifold limits have 
(relatively) few complex structure parameters. 

Another large class of examples can be obtained as P 1 bundles over toric surfaces: 

P 1 



B'. 

Let us start with the case B' = P 2 . The P^bundle over B' will be specified by an integer 
n according to the following toric data: 

D x D 2 D 3 D 4 D 5 
<D* 1 1 1 n (5.1) 
C* 1 1. 

D$ are the toric divisors : Di, -D3 are the pullbacks of the three lines in P 2 
and D4 and -D5 are the two sections of the P 1 fibration. The arithmetic genus of the 
fourfold divisors Vi = ir*(Di) is given by x(Vi, OyJ = l/2KsD^i = 1, . . . , 5, where Kb 
is the canonical divisor of B ||57|] . An immediate computation gives 

x(Vi) = x(v 2 ) = x(v 3 ) = -1, x(v 4 ) = -^f^, x(v 5 ) = n( ~ n 2 +3) . (5.2) 

Now we see that that the pullbacks to the fourfold of the two sections can not have 
simultaneously arithmetic genus 1. However, let us first choose n to be 1 or —1 so that 
either V4 or V5 will have arithmetic genus 1. Now we can enforce an ADE type singularity 
along the section whose pullback does not contribute initially to the superpotential. The 
components of the singular fiber will have arithmetic genus 1 and will project to that 
section. Since the two sections are linearly independent, they generate the Picard group 
of B and we have thus obtained a model where the condition for stabilizing all the Kahler 
moduli is satisfied. 

Note that in this case it is possible to stabilize the Kahler moduli using only divisors 
of gauge type, by enforcing ADE singularities of the Weierstrass model along both of the 
two sections of the P 1 bundle. However, this will not be true for models with h 1 ' 1 > 2. 
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We can perform a similar analysis in the case when B' is a Hirzebruch surface F m ; 
The toric threefold B is specified by two positive integers n and p as B = P(0F m © 
O-F m (nC +pf)), where C and / are respectively the negative section and fiber of F m . 
The toric data of B is 

£>i D 2 D 3 D 4 D 5 D 6 
C* 1 1 m p 
<D* 1 1 n 
C* 1 1. 

Di, % = 1, . . . , 6 denote again the toric divisors. Another quick computation gives 



(5.3) 



x(Vi) = 0, x(V 2 ) = 0, x(V 3 ) = -m, x(V 4 )=m, 

mn(n + l) mn(-n+l) ( 5 - 4 ) 
X(V 5 ) = -n-p-npH , x(V 6 )=™ + p-np . 

In order to satisfy the requirement stated above we need to have m = 1. It is easy to 
see that we can not have x(Vs) = xO^g) = 1- These equations would imply that n and p 
should satisfy (n —p) 2 —p 2 = 2 which however does not admit integer solutions. But since 
x(y$) an d x(Ve) are the the pullbacks of the two sections we can enforce again an ADE 
singularity of the Weierstrass model along one of them while choosing either n = 2(p + 1) 
or n = 2{p — 1) such that the pullback of the other section has also arithmetic genus 1. 
Obviously, in this case we can not construct a basis of the Picard group of B consisting 
only of arithmetic genus one divisors of gauge type, because h 1,1 = 3. 

A similar analysis can be carried in the case when B is a P 1 bundle over a toric 
del Pezzo surface. By making specific choices for the data of the P 1 bundle, in the case 
of the del Pezzo surfaces dP 2l dPs and dP^ we can construct examples that have a full 
basis of divisors of instanton type contributing to the superpotential. What happens if 
we also consider divisors of gauge type? For any P 1 bundle over dP 2 or <iP 3 , enforcing 
for example an ADE singularity along one of the sections, we obtain a model satisfying 
the above criterion. In the case of dP^ we have to enforce singularities of the Weierstrass 
model along both sections of the P 1 bundle. 

To summarize, there are several models with toric Fano base in which instanton-type 
divisors can stabilize Kahler moduli. In these models, which can be analyzed using existing 
techniques, the presence of a suitable nonperturbative superpotential is clear. 

There are also several possibilities for models with P 1 -fibered base in which gauge-type 
divisors can stabilize Kahler moduli. These models have heterotic duals and are potentially 
simpler, but establishing the existence of a suitable superpotential in these models requires 
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controlling the matter content and matter superpotential of the gauge theories. This is a 
rather complicated problem in the F theory framework, which has not been solved in the 
detail we need; in particular the flux contributions to the matter superpotential are not 
yet well understood. Thus, we will not reach definite conclusions for these models in this 
work. 

Finally, there are surely many more models in which the base is not Fano (any model 
with h ' > 10 is necessarily of this type), and there may also be models whose base has 
other fibration structures. 



5.1. The T\% model 



One of the examples from |57j is B = JF 18 , a toric Fano threefold |79|j8^ , p4]| . By our 
previous discussion, Z will be a double cover of B, branched along its canonical divisor. 
Thus Z can be realized as a quadric in Y = W(0^ 1S © Oj^ ia (K^ 18 )) . 

Note that Y is not a weighted projective space and its toric data is given by 
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(5.5) 



The generic divisor in the linear system | — Ky \ is the smooth Calabi-Yau variety Z 
and is defined by the equation 



W 2 + WU fa,...asll X i i + U2 E 9 bl bsU X 



bi 



0, 



(5.6) 



i=l 



(ai,...,a 8 )' i=l (bi,...,b&)> 

where the sums are taken over sets of positive integers a±, . . . , as,b%, . . . ,bs which satisfy 



-bi +b 7 + b 8 


= 2(- 


-ai + a-? + a 8 ) 


= 2, 


h 


+ b 2 


-b 4 


= 2(ai + a 2 


- 04) 


= 2, 


-bi + 64 + b 6 


= 2(- 


-ai + CI4 + as) 


= 2, 


bi 


+ b 3 


-b 6 


= 2(ai + a 3 


- a 6 ) 


= 2, 


-b 2 + 64 + ^5 


= 2(- 


-a2 + CI4 + as) 


= 2, 


b 2 


+ b 3 


-h 


= 2(a 2 + a 3 


- a 5 ) 


= 2, 


-63 + 65+^6 


= 2(- 


-a 3 + a 5 + a 6 ) 


= 2. 















(5.7) 
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Using the invariance under reparametrizations, we can set / ai ag = and the hypersur- 
face equation becomes 

8 

W 2 + U 2 9 bl ...bsi[ X ^ =°- ( 5 - 8 ) 

(6i,...,b 8 )' *=1 

The holomorphic involution O is simply W — > — W. All the third cohomology of Z is 
odd under the (pull-back) of and in particular 0*0 = — O, where O is the holomorphic 
three-form on Z. The third cohomology of the quotient B = Z/Cl is trivial since B is 
toric. Conversely, H 2 (Z,1L) and H 4 (Z,T) are even under O*. Therefore, all the complex 
structure and Kahler deformations remain in the spectrum []15| . 



Thus, we obtain a IIB orientifold compactification on B of the type we want. 
Topological analysis 

This can be done using standard toric techniques |J. The toric data of JF 18 is 



D x D 2 D 3 D A A, D 6 D 7 D 8 

C* -1 00000 1 1 

<D* 1 1 -1 0000 

C* 1 1 -1 , v 

C* 1 1 -1 ^ J 

C* 0-1 01 1000 

C* -10010100 

C* -1 1 1 0. 

We choose .Dj,i = 1,...,5 as basis for H 2 (B,~K). We have the linear equivalence 
relations D 6 = D 2 -D 3 + D 4 , D 7 = D 8 = ~D X + D 3 - D 4 + D 5 . 

The Kahler cone is not simplicial: although it is five dimensional, it has six generators 
which are given by 

fti=D 7 , n 2 = D!+D 6 + D 7l K s = D 1 +D 3 + D 6 + D 7 , TZ 4 = D 3 + D 6l 

(5.10) 

TZ 5 = D x - D 2 + D 3 + D 6 + D 7 , Tie = D x - D 2 + D 3 + 2D 6 + D 7 . 
In order to define a large radius limit and Kahler coordinates, we need to choose a simplicial 



decomposition of the Kahler cone and pick one of the subcones [24J]. We note that there 
exists a simplicial decomposition such that one of the subcones is generated by IZj, j = 
1, . . . , 5 and therefore we take the Kahler form to be 

5 

J = J2 t J U J = (*2 + h + h)D x - t 5 D 2 + (t 3 + U + t 5 )D 3 + (t 2 + t 3 + t 4 + t 5 )D 6 
j=i 

+ (t 1 + t 2 + t 3 + t 5 )D 7 . 

(5.11) 
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The dual polyhedron V z which encodes the divisors of the Calabi-Yau threefold Z 
has the following vertices: 

(0,0,-1,1), (0,1,1,1), (0,-1,0,1), (0,1,0,1), (0,0,1,1), 

(5.12) 

(0,-1,-1,1), (1,0,0,1), (-1,0,-1,1), (0,0,0,-1). 

The Hodge numbers! of Z are h l > l (Z) = 5, h 2 ^(Z) = 89. 

Next, we need to know how divisors of B pull back to the Calabi-Yau fourfold X. 
Since X elliptically fibered over the base B = JF 18 , we can construct it as a hypersurface 
in a toric variety. 

Using this, we can construct the dual polyhedron Vx, which encodes the divisors of 
the Calabi-Yau fourfold X. It has vertices: 

(0,0,-1,2,3), (0,1,1,2,3), (0,-1,0,2,3), (0,1,0,2,3), (0,0,1,2,3), 

(0,-1,-1,2,3), (1,0,0,2,3), (-1,0,-1,2,3), (0,0,0,-1,0), (0,0,0,0,-1). 

(5.13) 

Standard toric methods give for the Hodge numbers of the fourfold h ' = 6, h 2,1 = 
0, /i 3 ' 1 = 2194, h 2 ' 2 = 8844, x = 13248. We note that it is possible to find a triangulation 
of Vx consistent with its elliptic fibration structure such that each of the top dimensional 
cones having unit volume, thus guaranteeing smoothness of the corresponding Calabi-Yau 
fourfold. 

Nonperturbative superpotential 

The fibration n : X — > B has a section a : B — > X , n o a = 1b- The following 
toric divisors have arithmetic genus one: Vi = (0,0,-1,2,3), V2 = (0,1,1,2,3), V3 = 
(0,-1,0,2,3), V 4 = (0,1,0,2,3), V 5 = (0,0,1,2,3), V 6 = (0,-1,-1,2,3) and S = 
(0, 0, 0, 2, 3). The first six divisors are vertical, Vi = 7r*(L>i), i = 1, . . . , 6, and therefore 
contribute to the F-theory superpotential, while E does not contribute, since is the section 
of the elliptic fibration, E = cr(B). Moreover, since Vj,j = 1, ... ,6 are vertices of Vx, 
h°' l (Vj) = 0, i = 1,2, 3, j = 1, ...,6. Now, since D 1 ,...,D 5 generate B 2 {B,TE), we see 
that this model has a complete set of contributing divisors and satisfies the requirement 
of section 3. 

Note that the divisor Di on the base B is the exceptional divisor that corresponds to 
the contractions of both first and sixth Mori cone generators. This is possible since D\ is 
isomorphic to a Hirzebruch surface Wq which is a product of those two curves. 

7 We have used the computer program POLYHEDRON, written by Philip Candelas. 
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In order to study the question of fixing the Kahler moduli, we need to compute the 
volumes of the divisors contributing to the superpotential, as well the volume of the three 
dimensional base B. To achieve that, we triangulate the fan of B. We list the cones below: 



£>i£> 4 £>8, DxDqDt, DiDeDs, D 2 D 4 D 7 , D 2 D 4 D 8 , 

(5.14) 

D 2 D 5 D 7 , D 2 D 5 D 8 , D 3 D 5 D 7 , D 3 D 5 D 8 , D 3 D 6 D 7 , D 3 D 6 D 8 . 

Using the above (unique) triangulation, we obtain the following nonvanishing inter- 
section numbers 



Df = 2, 


D 2 D 4 = -1, 


D 2 D 6 = -1, 


DjD 7 = -1, 


D 3 - 
u 2 — 


-1, 


D 2 2 D 4 = 1, 


D\D 7 = -1, 


Dl = -1, 


£>lA> = 1, 


D 2 D 7 


= -1, 


D 2 D 2 = -1, 




D\D 7 = -1, 


D 2 D 2 = 1, 


DlD 3 


= 1, 


Dl = -2, 


D 2 5 D 7 = -1, 


D 2 6 D 3 = -1, 


= 1, 


D 2 D 7 


= -1. 












(5.15) 


We obtain for the divisor volumes Tj 


= DiJ 2 /2 and the total volume V 


= J 3 /6: 





Ti = t!U, r 2 = -|(2t! + 2t 2 + 2t 3 + t 5 ), 

T 3 = |(2t! +t 2 + 2t 3 + 2t 5 ), T 4 = ^(t 2 + t 3 )(2ti + *2 + * 3 ) 

T 5 = (*3 + t 4 )(tl +t 2 +t 3 + t 5 ), Tq = i(t 3 + *5)(2*1 + *3 + *5), 

tota t\tk 9 to to^4 tq^4 ^ 

U = ti* 2 *3 + + + *2*§ + 4 + *1*2*4 + + *1*3*4 + *2*3*4 + ^ 

t 2 t§ t 2 t 2 

+ tit 2 ts H — - — h tit 3 t 5 + 2t 2 t 3 t$ + t 3 t 5 + tit 4 t$ + t 2 t 4 t$ + t 3 t 4 t$ H — — 

2 2' 

Complex Structure Moduli 

To compute or count flux vacua in arbitrary regions of the complex structure moduli 
space of the Calabi-Yau threefold Z, one would have to compute the periods for the 
generic hypersurface Z. These are generalized hypergeometric functions in 89 variables. 
In principle they can be computed using existing techniques, but this would require a lot 
of work, even using a computer. 

A somewhat easier (but still formidable) task is to describe the periods in the vicinity 
of the large complex structure limit. This is equivalent to computing the triple intersections 
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for the mirror threefold Z, which has h 2,1 = 5 and h ' = 89. To make the description of 
this part of the moduli space complete, one furthermore has to compute the Kahler cone, 
i.e. the part of the parameter space IR 89 where all holomorphic curves have positive area. 
We describe the algorithms we used to achieve this in appendix B. 

Flux vacua 

Counting flux vacua can be done using the techniques of p| j29|1 as outlined in section 
4. Since L* = 13248/24 = 552 ^> 63 = 180, we expect that the approximations made to 
derive the counting formula (|4.1|) should be valid. 



According to this formula, the expected number of vacua with cosmological constant 
e K °\W \ 2 less than A* < L*, equals (2ttL) &3 / & 3 ! ~ !0 307 multiplied by K/L* times the 
integral of a geometrical density function. Taking the integration domain M. c equal to the 
entire moduli space, the estimate ( |4.2|) indicates that the geometrical factor should be of 
the order of the Euler characteristic of Ai c . This does not need to be an integer, as the 



moduli space is a noncompact orbifold, and in examples [|29[ is a small fraction, of order 
l/|r|, where V is the order of a finite group (e.g. 2 5 for the mirror quintic) or the volume 
of a group [771. This will be far larger than 10 -307 and thus we expect many vacua with 



very small cosmological constant. 

Can we find such vacua at large complex structure, or equivalently, at large volume 
of the mirror? We can judge this by integrating the density function over, say, the region 
defined by requiring all curve areas of the mirror to be bigger than 1. Using the approach 
for estimating the density function po explained in appendix A and the construction of the 
geometry of the moduli space outlined in appendix B, we have done Monte Carlo estimates 
of this integral. The results are as follows. The average value of \i defined in appendix A 
is of order 10 90 , hence the density function po is of order 10 240 , and we take this out of the 
integral ( }4.1| ). The remaining volume integral, evaluated using 10 7 Monte Carlo sample 
points, gives a number of order 10~ 650 . Putting everything together we get, up to the 
factor A*/L*, 

N vac (LCS) ~ 10" 100 . (5.17) 

Why is the volume so small? One important reason is the mirror volume suppression 
factor V~ n / 3 mentioned at the end of appendix A. For curve areas y l bigger than one, 
the mirror volume V > 10 10 . This becomes understandable when one considers that when 

V is written in terms of the curve area coordinates y l (cf. appendix A), the expression 

V = D^-\y l yiy k /6 contains ~ 10 6 terms, and the D^jl are widely distributed between 
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and 10 8 . With V > 10 10 , the volume suppression factor is ~ 10 300 . Additional 
suppression comes from the fact that det g ~ 10 -100 for typical values of y with V(y) = 1, 
and from the smallness of the Euclidean volume of the surface V(y) = 1. It is possible 
that the Monte Carlo missed regions in which V is smaller than 10 10 , or that some of this 
is an artifact of our approximate parametrization of the Kahler cone with the curve areas 
y l , but we see no particular evidence for this. 

Thus it seems likely that any flux vacua in this region are special cases, and there is 
no reason to expect a multiplicity of vacua out of which some would have small Wq. This 
paucity of vacua in the large complex structure limit is not specific to this example, but 
rather is a very general feature of models with many moduli, as explained in ||30|| . 

If we neglect world-sheet instanton corrections, then by lowering the cutoff to y l > 
0.075, the expected number of vacua becomes of order 1. In few modulus examples, the 
instanton sums tend to converge all the way down to zero, so this might be a valid indication 
of where vacua will start to exist. 



5.2. The T\\ model. 



Another example from Grassi's list is the base B — T\\ [83,83] . We now construct an 
explicit toric model of a Calabi-Yau fourfold elliptically fibered over this base. This model 
has an orientifold limit Z with h ' = 3, h 2,1 = 111. 

The toric data of T\\ is given by fF3|,[F9],|83], 





D t 


D 2 


#3 


D 4 


D 5 


D 6 


<D* 





-2 


1 


1 


1 





<D* 


-1 





1 








1 


<D* 


1 


1 


-1 








0. 



(5.18) 



The generators of the Kahler cone of B are 



TZ t = D 5 , TZ 2 



D 6 , 



1lz=D l + D 6 . 



(5.19) 



We have the linear equivalence relations D 4 = D 5 = Di + D 3 and D 6 = Di+D 2 +2D 3 . 
The Calabi-Yau fourfold elliptically fibered over B may be constructed as a hypersurface 
in a toric variety. The dual polyhedron V, which encode the divisors of the Calabi-Yau 
fourfold has vertices: 

(1,1,0,2,3), (0,-1,0,2,3), (1,0,0,2,3), (0,-1,1,2,3), (-1,-1,-1,2,3), 
(0,1,0,2,3), (0,0,0,-1,0), (0,0,0,0,-1). 
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The Hodge numbers of this fourfold are h 1 ' 1 = 4, /i 2 - 1 = 0, /i 3 ' 1 = 3036, /i 2 ' 2 = 12204, x = 
18288. Again, we find that the fourfold is smooth. 

The fibration n : X — > B has a section a : B — > X , n o a = 1b- The following divi- 
sors have arithmetic genus one: Vi = (1, 1, 0, 2, 3), V 2 = (0, -1, 0, 2, 3), V 3 = (1, 0, 0, 2, 3) 
and S = (0,0,0,2,3). The first three divisors are vertical, Vi = 7t*(Di), i = 1,2,3, and 
therefore may contribute to the F-theory superpotential, while £ does not contribute, 
since is the section of the elliptic fibration, £ = o~(B). It is possible to check that in 
fact h°' l (Vj) = 0,i,j = 1,2,3, thus all the vertical divisors do give a contribution to the 
superpotential. Therefore, this model also provides a complete basis of divisors. 

To compute the volumes of the divisors Di, D2, -D3, we triangulate the fan of B. 
Its 3-dimensional cones are DiD^Dq, DiD 5 D g , D^D^Dq, D 2 D^D 5 , DiD^D^, DiD 3 D 5 , 
D2-D3-D4, D2-D3-D5. Using this (unique) triangulation, we obtain the following nonvanish- 
ing triple intersections: 

D 3 = -3, D 2 D 3 = 2, D 1 D 2 3 = -1 1 D| = 4, D 2 2 D 3 = -2, D 2 D\ = 1. (5.20) 

Let J = X^i=i tiTZi be the Kahler form of B. We obtain for the divisor volumes 
T t = DiJ 2 /2 and the total volume V = J 3 /6: 

t t 2 

7-i = -|(2ti+t 2 +4t 3 ), T 2 =-±, T 3 = t 3 (t 1 +t 3 ), 

(5-21) 

y = tjh + tlt 1 + t^ + t 1 t 3+ + t a tg + ^2 + 2hf 2 + _*3 

5.5. T/ie orientifold o/P [1,1,1,6,9] 

Finally, out of the various possibilities which can be obtained as P 1 bundles over 
a toric surface, the simplest is perhaps the P 1 bundle over P 2 . The toric data of the 
threefold is presented in ( |5.1| ), where we take n = —6, so that B = P((9 P 2 © 0-pz{— 6)). 
The toric data for B is as follows: 

D x D 2 D 3 D 4 D 5 
C* 1 1 1 -6 (5.22) 
<D* 1 1. 

The toric divisors D4 and -D5 are the sections of the P 1 bundle over B' = P 2 . We 
will construct an elliptic fourfold X over B in such a way that X will be the resolution 
of a Weierstrass model W which has a D4 singularity along the first section and an Eq 
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singularity along the second section. The motivation for this choice will become clear later, 
when we study Kahler moduli stabilization for these models 

The dual polyhedron Vx, which encode the divisors of the Calabi -Yau fourfold has 
vertices: 

(-1,-1,6,2,3), (0,0,-3,2,3), (0,1,0,2,3), (1,0,0,2,3), (0,0,-2,1,1), 
(0,0,-1,0,0), (0,0,0,-1,0), (0,0,0,0,-1). 

The Hodge numbers of this fourfold are h 1 ' 1 = 13, h 2 ' 1 = 0, h 3 ' 1 = 1071, h 2 ' 2 = 4380, x = 
6552. Again, we find that the fourfold is smooth. We can also enforce an E? or E% 
singularity along the infinity section. The data for the fourfolds obtained this way are, 
respectively, h 1 ' 1 = 14, h 2 * 1 = 0, h 3 ' 1 = 935, h 2 > 2 = 3840, \ = 5742 and h 1 ' 1 = 16, h 2 ' 1 = 
253, /i 3 ' 1 = 745, h 2 > 2 = 2582, X = 3096. 

It is easy to check that in this case the orientifold limit Z will be an elliptic fibration 
over P 2 , which is familiar as the hypersurface in weighted projective space Pm 1169 j 



studied in [|21| , |33H and several other works. In order to do that, we note that Z is given by 



a quadric in a toric variety Y described by the following data 

X 1 x 2 x 3 x 4 x 5 u w 

C* 1 1 1 -6 3 , v 

C* 1 1 -2 [ } 

C* 1 1. 

This threefold has h 1,1 = 2 and h 2,1 = 272. There is a unique toric Calabi- Yau threefold 
with this Hodge numbers, the elliptic fibration over P 2 17ej 



Nonperturbative superpotential 

We proceed as before and start by listing all toric divisors of arithmetic genus one for 
the first of the models above. These are: Vi = (0, 0, -3, 2, 3), V 2 = (0, 0, -2, 1, 1), V 3 = 
(0,0,-2,1,2), V 4 = (0,0,-2,2,3), V 5 = (0,0,-1,0,0), V 6 = (0,0,-1,0,1), V 7 = 
(0, 0, -1, 2, 3), V 8 = (0, 0, 2, 2, 3), V 9 = (0, 0, 1, 2, 3). 

The first seven divisors project to D$ and the last two project to D^. These are going 
to be the only base divisors contributing to the superpotential for the other two models as 
well. 

To compute the volumes of the divisors -D4 and -D5 we start by triangulating the fan of 
B. Its three dimensional cones are D1D2-D4, -D1-D2-D5, -D1-D3-D4, D1D3D5, D2-D3-D4 and 
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L>2 -D3-D5. Using this triangulation we obtain the following nonzero intersection numbers 

DjD 4 = l, DjD 5 = l, D 2 D 4 = 1, D 2 D 5 = l, D 2 D 4 = 1, D 2 D 5 = 1, 

DiDj = -6, DxD 2 5 = 6, D 2 D| = -6, D 2 D 2 5 =6, D 3 D 2 4 = -6, D 3 D 2 = 6, (5.24) 

D| = 36, D5 = 36. 

Let J = ti-Di + t^Dc, be the Kahler form of B. We obtain for the divisor volumes 
n = DiJ 2 /2 and the total volume V = J 3 /6: 

*? (ti+6t 5 ) 2 
T4 = — , r 5 = , 

2 2 (5.25) 

V= -(3tJt 5 + 18Ml+36ti). 
o 

Complex structure moduli and flux vacua 

The full 272 parameter prepotential for this model has not been worked out. However, 
this moduli space admits a T = Zg x Zis action, which fixes the two parameter subspace 
of CY's with defining equation 

/ = x\ 8 + xl, 8 + x\ 8 + x\ + x\ — \%i\)X\X2Xj,x^x^ — ?><\)x\x\x\. 



This is the subset of CY's obtained by the mirror construction of [|58| and the six periods 



of these mirror CY's, which is the same as the subset of periods of T-invariant cycles, are 



worked out in [21)] . 



If one turns on flux only on these cycles, since the resulting superpotential and Kahler 
potential are T-invariant, one is guaranteed that DiW = in all T-noninvariant directions, 
and thus one can find flux vacua just by working in the T-invariant part of the moduli space, 
call this M.Y- (This observation is also made in fl53|.) Of course, one would eventually 
need to check that all moduli remain non-tachyonic after supersymmetry breaking. 

For the choices of gauge groups we considered, L = x/24 ~ 100 — 300, so the total 
index for flux vacua in this subspace is P,p9 



(2ttL 



16 



6! 



/ det(-R-u). 

JMrxn 

As mentioned earlier, the integral can be estimated as 1/12 |T| = 1/1296, leading to 

7~ 10 11 - 10 13 . 
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(5.26) 



The number at weak coupling 1/g 2 < e will be roughly e times this, and the number with 
e K °|Wo| 2 < ALT 3 will be roughly A times this, so there are clearly many weakly coupled 
flux vacua with small Wq. 

Just to get a few explicit flux vacua, we consider the region of large complex structure, 
i.e. the region in which instanton corrections are small, Nie 27riWi « 1. From [El, the 
prepotential at third order in the world-sheet instanton expansion is 

T =- (9to? + 9w\w 2 + Sw^l) 
6 

9 2 3 17 3 

- - 2 WlW2 ~ ~4 Wl ~ 2 W2 + * 

+ 7^3 ( 54091 + 3q2 + " 1080^2 - ^-q\ 

{2m) A \ 2 2 

+ 560qf + 143370^92 + 2700qiq| + ^4 + . . 

9 

with q, = e 2m and f = C 2 ( §ffl w -1.30843 i. 

Looking at the coefficients in this expansion, the instanton corrections should be small 
for W2 » 1/6 and w\ > 1. To find quantized flux vacua, we used the procedure discussed 
in section 4. As an example, we first look for fluxes which stabilize the moduli at the 
rational point r = 3z, w\ = i,W2 = i, ignoring world-sheet instanton corrections and 
rationally approximating £ by —13/10. One finds a quantized flux vacuum with W = 0, 
with fluxes (Nrr; Nns) equal to 

{0, 69, 28, 0, 0, -20; -49, -18, -6, -4, 0, 0}; L = 352 

granting that the quantization condition on the orientifold is the same as on the original 



CY. This was argued in general in [[47]]; the fact that some cycles are smaller on the 
orientifold which naively doubles the Dirac quantum, is compensated for by the possibility 
of discrete RR and NS flux at the orientifold fixed points. A careful discussion in the 



present example might be possible using K theory |8(| . 

Restoring the exact value of £, including the first instanton corrections and solving 
the resulting equations DW = produces a W 7^ vacuum near the starting point. We 
find 



r 




2.945z 


Wi 




0.9625i 


w 2 




1.1037i 



*K\W\ 2 = 1.379 x 10~ 4 
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Examples with larger r, Wi can easily be found in this way, but L tends to become bigger 
then as well. 

For these vacua, one can check that the instanton corrections are small (as one would 
guess by the small corrections they lead to for the moduli) and the vacua appear sound, 
on the level we are working. 

6. Numerical results on Kahler stabilization 

6.1. JF 18 model 

The Kahler moduli are stabilized for generic order 1 values of the bi. Taking Wq = 
10" 30 , the typical values for the U are £i ~ 50 - 100, U ~ 0.1 - 0.3 for i = 2,4, 5 and 
£3 ~ 0. The corresponding values of the contributing divisor volumes are Vi ~ 11 — 12, 
and the total volume V ~ 5 — 10. Taking the bi all equal gives a nongeneric singular 
solution: t\ = 00, ti = for i > 0. But with bi = (1, 1.5, 2, 2.5, 3, 3.5) for example, we get 
U = (53.3,0.209,0.00156,0.208,0.222), n = (11.1, 11.9, 11.2 + ztt, 11.2, 11.3 + ztt, 11.9 + m), 
and V = 7.32. Different choices for the phases of Wq and five of the bi can be absorbed in 
shifts of the imaginary parts of the (the axions). Note that since the areas of the Mori 
cone generators are (£1, £2 + £3, £4 + £3, £5 + £3, £5, £4, £2), these solutions lie well inside the 
Kahler cone; £3 = is only an interior boundary of a subsimplex of the Kahler cone. 

The typical values for the ti and the total volume are rather small for this model, 
so a' corrections could become important. (The situation is better for the other models 
we consider.) However, since the solution in terms of the Ti is mainly determined by the 
exponential factors (hence the near-equal values of Ti ~ — (In Wq)/2tt 11), it is reasonable 
to believe that it will still exist at approximately the same values of even after taking 
into account such corrections. 

Finally, to have a stable vacuum after lifting the cosmological constant to a positive 
value, the critical point has to be a minimum of the potential. We verified numerically 
that this is indeed the case. 

6.2. T\\ model 

The Kahler moduli are stabilized for generic order 1 values of the bi. For Wq = — 10~ 30 , 
bi = 1, we get ti = (4.89, 1.30, 1.76). These are also the areas of the three generators of 
the Mori cone. The corresponding complexified volumes of the divisors Di contributing to 
the superpotential are Tj = (11.8, 11.9, 11.7) and the total volume is V = 93.3. The critical 
point is a minimum of the potential. 
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6 - 3 - F fl, 1,1,6,9] modd 

Assuming the gauge theory generates a superpotential W = Wo + Yli=i bie~ 27zaiTi , 
where t±, t 2 are the complexified volumes of the divisors D4 and -D5, we find that the 
Kahler moduli are stabilized for generic order 1 values of the 6j, provided a± > a<i- The 
following are some of the values of ti, Vi and the volume V, obtained for bi = 1 and 
different choices of a%, a 2 and Wq: 

ai a 2 W h t 5 V 1 V 2 V 

1/4 1/30 1(T 30 9.83 2.76 48.3 348 484 . , 

1/4 1/30 10" 5 4.61 1.14 10.6 65.7 39.1 ^ ' 

1/4 1/12 10" 30 9.73 1.16 47.4 1 39 1 03 

1/4 1/12 10" 5 4.40 0.468 9.64 25.9 8.01 

The chosen values of ai correspond to pure G2 x E% resp. G<i x Eq gauge theory. Again, 
the critical point is a minimum of the potential. 

The fact that a± > 0,2 is needed to have a solution (lying well inside the Kahler cone) 
can be seen directly from ( |5.25| ): the approximate solution is T\ = — In Wo/2nai, but ( |5.25| ) 
implies T\ < T2, so a\ > a<i- 



7. Conclusions 

Our primary result is that we have candidate lb orientifold compactifications in which 
nonperturbative effects will stabilize all complex structure, Kahler and dilaton moduli. 

Models which stabilize all Kahler moduli by D3-instanton effects are not generic, but 
not uncommon either. We listed all the possibilities with Fano threefold base, and many 
possibilities whose bases are P 1 fibrations. It turns out that these models must have 
several Kahler moduli and several non-perturbative contributions to the superpotential, as 
in the early racetrack scenarios for moduli stabilization. 

We see no obstacle to adding supersymmetry breaking effects such as the antibrane 
suggested by KKLT, D term effects or others. One also expects these potentials to have 
many F breaking minima (statistical arguments for this are given in |29|). On general 



grounds, since the configuration spaces parameterized by other moduli such as brane and 
bundle moduli are compact, after supersymmetry breaking all moduli should be stabilized. 

The models with Fano threefold base are rather complicated, with many complex 
structure moduli. It is not clear to us that this makes them less likely candidates to 
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describe real world physics, but it is certainly a problem when using them as illustrative 
examples. 

There may well be simpler models among the P 1 fibered models, whose nonperturba- 
tive effects have a simple gauge theory picture. Perhaps the simplest is the P 1 fibration 
over P 2 or "11169 model." While the standard F theory analysis of the models we dis- 
cussed have suggests that they have too much massless matter in one gauge factor to 
produce non-perturbative superpotentials, we suspect that other effects, in particular flux 
couplings to these brane world- volumes, would lift this matter, leading to working models, 
and intend to return to this in future work. 

Another possibility for finding simpler models would be to stabilize some of the Kahler 
moduli with D terms. Following lines discussed in ||, one can show that in configurations 
containing branes wrapping k distinct cycles, D terms will generally stabilize k — 1 relative 
size moduli. Since this relies on a' corrections, the resulting values of Kahler moduli will 
be string scale, but order one factors in the volumes and gauge field strengths might be 
arranged to make this a few times the string scale, which could suffice. 

In any case, even the simplest models under discussion have many more moduli to 
stabilize. Furthermore, one must check that a candidate vacuum is not just a solution, 
but has no tachyonic instability. Now once one has established that these moduli indeed 
parameterize a compact configuration space, it is clear that the minimum of the effective 
potential on this space will be a stable vacuum. While it may be hard to compute the value 
of the potential at this minimum, the large number of flux vacua strongly suggests (as in 
||13|| ) that whatever it is, it can be offset by a flux contribution to lead to a metastable de 
Sitter vacuum. Thus, granting the effective potential framework, it will be extremely sur- 
prising if vacua of this type do not exist, while technically quite difficult to find them. This 



is the type of picture which motivates the statistical approach discussed in |4^,p|,|7|j29|,p0 
as well as anthropic considerations such as []8|,|9T|,f7i ■ 

Can we say anything about the validity of effective theory? At present we see no 
clear reason from string/M theory or quantum gravity to doubt it. However, even within 
the effective field theory framework, there is another important assumption in KKLT, our 
work, and the other works along these lines. Namely, we have done Kaluza-Klein reduction 
in deriving the configuration space of Calabi-Yau metrics, and typically will take similar 
steps throughout the derivation of the effective field theory. Could it be that in many of 
these backgrounds, some of the KK and stringy modes which are dropped in this analysis 
are in fact tachyonic? Since these constructions rely on approximate cancellations between 
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many diverse contributions to the vacuum energy, it is conceivable that subsectors of the 
theory have instabilities which do not show up in the final effective lagrangian; this should 
be examined. In any case, a lot of work remains to see whether these models are as plentiful 
as they now appear. 
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Appendix A. Counting flux vacua when 63 is large 

Equation (|4.1|) gives the expected number of flux vacua with L < L*, \Wq\ 2 < A 



t e H C H and z G M. c C A4 C , for A* <C L*. Here we will derive expressions for p 
suitable for Monte Carlo estimates in the case of many moduli. 
It is somewhat more convenient to rewrite (14.11) first as 



N vac = vol(H) / d An z det<7 u(z). (A.l) 

>M C 



The density function u(z) is |29 



u(z) = ^ [ d 2 X [ d 2n Z I det ( I ^% fc ) | 2 + 0{X). (A.2) 



The integration variables X and Z l range over C and are the remnants of the flux vector 
after a certain ^-dependent change of basis and after imposing the constraint DWq = 
(now considered as a constraint on the continuous flux components at a given point z). 
Indices are lowered with the metric g^. The factor in front of the integral comes from 
the Jacobian corresponding to the change of basis, and the determinant in the integrand 
is the I det .D 2 Wo I Jacobian accompanying the delta function imposing DWq = 0. Finally, 
the Tijk are the "Yukawa couplings" characterizing the special geometry of the complex 
structure moduli space, i.e. 

Fijk = e Kc / A didjdkQ. (A.3) 
Jz 

37 



Since we are interested in very small values of the cosmological constant, i.e. A* <C 1, the 
0(X) part can be dropped from the integrand and the integral over X simply gives a 
factor 7rA*. 

By rewriting the determinant in the integrand as a Gaussian fermionic integral and 



the bosonic integral as a Laplace transformed Gaussian ||29| , and then doing the bosonic 
integral, this expression can be rewritten as a fermionic integral with quartic fermionic 
action, hence it reduces to a finite number of terms. (Alternatively, one could stick to 
bosonic variables and apply Wick's theorem.) This is useful for low n, but for large n the 
number of terms becomes enormous, about n 4n , and straightforward numerical evaluation 
becomes impossible. Instead, we want to rewrite ( |A.2| ) in a form suitable for Monte Carlo 
estimates. 

To this end, we implement the constraint \\Z\\ 2 < L* by inserting f^* d£5(£ — \\Z\\ 2 ) 
in the integral. Changing variables from Z to U with Z l = ^flU 1 , and doing the ^-integral 
then gives 

9 2rc+2 \ r2n+l /• / n tj \ 

"<'> = (2„ + y d ; tg /^OT*-l)|de t (° JfrF- (A ' 4> 
Now define the following "spherical" average, for any function f(U): 

(f(U)) Jd^US(\\Uf-l)f(U) 



Noting that 



we can thus write 



with 



d 2n U 8(\\U\\ 2 - 1) = - -— (A.6) 

v " " ' (n-l)!det# v ' 

o2n+2_n+l\ r2n+l 



11^11=1 



Comparing this to (|4.1|) , we get 



, N (2n + 2)(2n)! , , . 

This can be made more explicit in the large complex structure limit of Z, parametrized 
with the special coordinates t % = x l + iy\ where y l becomes large and Xi G [0, 1]. All data 
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is encoded in the triple intersection numbers Dijk of the mirror Z to Z. The Kahler 
potential is KL C = — ln(2 3 F), where V = ^Dijky l y-'y k ■ The metric and its determinant are 



ViVi Va , (-l) n+1 detl/ i? 

Hj 41/2 4 y' y 2 2n+1 l/ 



det ^= L,- (A.10) 



where Vi = diV = ^Dijky^y h and Vij = didjV = Dijky h . Finally, the Yukawa couplings 
are Tijk = ^"Dijk- Using all this and pulling the factor out of the determinant, (|A.8| ) 
becomes 

Note that n(y) is invariant under rescaling y l \y l . 

The average (f(U)) can be estimated numerically using Monte Carlo methods; for 
many variables this is in fact the only possible way. In the simplest version, one repeatedly 
picks a random vector Z from a normal distribution with mean zero and covariance matrix 
proportional to the metric gq (or any other distribution depending only on ||Z||), one 
evaluates f(U) with U = Z/\\Z\\, and at the end one computes the average of the values 
obtained. This gives an approximate value for po(y), the approximation becoming better 
with increasing number of sampling points. 

What remains to be computed then to get the number of flux vacua is the integral 
Again this can be done by Monte Carlo integration. However, computing po(y) at 
every sample point in the way described above would require far too much computation 
time in the models we consider. Therefore we will replace po(y) by its average over a 
limited number of sample points. This is a reasonable thing to do, since in examples po(y) 
stays of more or less the same order of magnitude over the large complex structure part 
of moduli space (in particular, as noted before, it is scale invariant). The integral then 
becomes proportional to the volume of the part of moduli space under consideration, which 
can estimated numerically in reasonable time. 

One important universal property of this integral can be deduced directly: vacua are 
strongly suppressed at large complex structure. To see this, note that under rescaling 
y — ► Ay, d n y det g — > X~ n d n y det g, so if for instance the region under consideration is 
given by V > V*, then 

N vac (V > V,) oc 1/V* n/3 . (A.12) 



We discuss this in more detail in 3C 
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Appendix B. Taming the complex structure moduli space of the JF 18 model 

The Calabi-Yau 3- fold Z of the model of section 5.1 has 89 complex structure 
moduli. Describing this space completely including its exact periods would be extremely 
complex, so we restrict ourselves to the large complex structure limit, which can be con- 
structed as the classical Kahler moduli space of the mirror Z. Even this poses quite a 
challenge. 

The classical Kahler moduli space of Z is specified entirely in terms of the triple 
intersection numbers of the divisors Di . These in turn can be computed from the quadruple 
intersection numbers of the divisors of the ambient toric variety Y . There are more than 
100 independent divisors in Y, so a priori there are more than 100 4 = 10 8 intersection 
numbers to compute. 

Before one can start doing this, one needs a maximal triangulation of the fan of the 
toric variety. Although the fan in this case is very large (116 points), we managed to do 
this by hand, and found 576 cones of volume 1. 

The next step is to determine the actual Kahler cone within the IR 89 parameter space, 
i.e. the cone in which all holomorphic curves have positive area. This involves constructing 
and solving 377 inequalities in 89 variables, which is more difficult than one might expect. 
Even finding just one point satisfying these inequalities takes several hours on a PC. 

In the end we want to compute numerically the volume of a subspace of the Kahler 
cone, as described in appendix A. Because of the high dimensionality of the space, this 
needs to be done by Monte Carlo. The integrand is relatively costly to evaluate, as it 
involves computing the 89 x 89 determinant given in (|A.10| ). On a 2.4 GHz Pentium, 
evaluation using 10 7 Monte Carlo sample points took about three days. 

In the following we outline how we proceeded to achieve these goals. 

B.l. Intersections 

The toric variety Y has 116 toric divisors Di, which correspond to points p l G IR 4 in 
the polyhedron V-p. The vertices of Vy are given by 



(4,0,-2,1), (-2,-2,0,1), 
(-2,0,2,1), (0,0,2,1), ( 
(-2,-2,2,1), (0,-2,2,1), 



(-2,2,0,1), (0,0,0,-1), (2,2,0,1), 
2,2,-2,1), (4,2,-2,1), (2,-2,0,1), 
(-2,0,-2,1). 



The triangulation of the polytope has 576 cones DiDjDkDi, all of volume 1. 



40 



Fig. 1: Partial triangulation of a 3- face of the polyhedron. Black dots correspond 
to vertices, blue dots correspond to points lying on edges or codimension two faces, 
while purple dots correspond to points interior to the 3-face. 

The rule for the intersection product of four distinct divisors is the following: if the 
four divisors span a cone, their intersection is 1, otherwise it is 0. It is more complicated 
to find intersections where some divisors are the same, such as ■ ■ Dj ■ D^. This is 
done by making use of the 4 linear equivalence relations that exist between the 116 toric 
divisors, which can be simply read off from the points: p\Dt = 0. This gives four equations 
p Dnjk = for each distinct triple ijk. The unknowns are the double index intersection 
numbers Dajk, so there are (more than) enough equations to find these. Solve this using 
a computer as one system of millions of equations in millions of variables would take a lot 
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of time and memory. However, the problem can be split up in a much smaller number of 
systems of four equations in three variables: for a given i < j < k, the three variables are 
Dujk, Dijjk, Dijkk, and the only ijk which need to be taken into consideration are those 
for which DiDjDk actually appears somewhere as a face in the list of cones. For other 
ijk the three unknowns are trivially zero. A similar reasoning can be followed to compute 
intersections with three and four identical indices. Thus almost all intersections vanish, 
and the remaining ones can be computed by computer in less than a minute. 

Let us now turn to the Calabi-Yau hypersurface Z in Y. Of the 116 divisors Di in Y, 
only 93 intersect Z and descend to divisors on the Calabi-Yau (89 of those are independent). 
We will denote these divisors by the same Di. Their triple intersection numbers are now 
easy to compute. They are simply given by the intersection with the anticanonical divisor: 
Dijk = Dijki- The resulting volume function is V = \Di^x l x^ x k . 

Finally, the large complex structure prepotential of the original Calabi-Yau Z is given 
by this expression with the x % replaced by complex variables, and takes the form 

T = t\ - t 2 2 1 3 + 780 more terms. 

It is available upon request. 

B.2. Mori cone 

A class in H 2 (Z) is specified by its intersections with the divisors Di. It is important 
to know which of these classes can effectively be realized as holomorphic curves. The set 
of such effective curves forms a cone, called the Mori cone. The generators of the Mori 
cone C , given by their intersection numbers C a i with the divisors Di, are contained in the 
list of 'special' linear relations C a ip % = between the polyhedron points p l corresponding 
to the divisors. There is one such special relation for each adjacent pair of cones in the 
triangulation of the polyhedron, found as follows. Denote the three common points of the 
two cones by /i, f 2 and fs, and the two additional points by p and q. Then there will be 
a relation p + q + n\f\ + n 2 f 2 + ^3/3 = 0, with the integer. 

Applying this to our model gives 579 relations. These correspond to effective curves 
of Y. What we want however are effective curves of Z ', and these are obtained by keeping 
only the curves having zero intersections with the 23 divisors of Y which do not descend 
to divisors of Z. This reduces the list to 377. Only a subset of this list will constitute a 
basis of generators of the Mori cone: any curve that is a positive linear combination of 
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the others can be dropped. A basis can thus be constructed in steps as follows. Start 
with the first curve. Add to this the second one and check if in the resulting set any one 
of the curves is a positive multiple of the other. If so, remove this curve. Then add the 
third curve and remove any curve in the resulting set that is a positive linear combination 
of the others. Then add the fourth curve, and so on, till all 377 candidates have been 
considered. This procedure is much faster than starting with all curves and removing the 
dependent ones one by one, because verifying if a given vector equals some positive linear 
combination of a set of n vectors becomes nontrivial if n is bigger than the dimension, and 
with n substantially bigger, it is computationally extremely expensive (hours for n — 377 
in this case). 

The resulting basis of the Mori cone consists of 111 generators C a . 
B.3. Kahler cone 

The dual of the Mori cone is the Kahler cone, that is the set of all J = x l Di such 
that C a ■ J > 0. Note that this is a tiny fraction of the M 89 parameter space; an estimate 
for the probability of a random point to be in the Kahler cone is (1/2) 111 ~ 10 -34 . Since 
we want to integrate over the Kahler cone, it is important to have a good parametrization 
of it — multiplying the integrand by a step function with support on the Kahler cone 
and Monte Carlo integrating over all x certainly won't do the job, since effectively all 
sample points will integrate to zero. Ideally, one would construct the exact generators 
K p of the Kahler cone and write J = t p K p , t p > 0. There exists an algorithm to find 



these generators, decribed in Jl8[ p.ll. Unfortunately, because the Mori cone is far from 
simplicial, this involves running over ( 8 g ) ~ 10 23 candidate generators (the rays obtained 
by intersecting 88 of the 111 zero planes). It would take about the age of the universe to 
complete this task on a PC. Nevertheless, with some luck, we were able to construct an 
approximate parametrization. Applying Mathematica's function Inequalitylnstance to 
find one solution to our set of 111 inequalities results in a point which has exactly 89 curve 
areas Cj ■ J equal to 1, and 22 bigger than 1. Taking the areas of these special curves as 
coordinates yi, it turns out that for uniformly random positive y values, the resulting J 
generically lies inside the Kahler cone. The integral over the Kahler cone can therefore be 
done by Monte Carlo integrating over all positive y and multiplying the integrand by a 
step function with support on the points satisfying the remaining 22 inequalities. In this 
way, not too many sample points evaluate to zero. 
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